arXiv:1509.05834v2 [math.OC] 26 Sep 2015 


Control of the Landau—Lifshitz Equation 


A.N. Chow and K.A. Morris 

Department of Applied Mathematics, University of Waterloo, Canada 
a29chow@U'waterloo. ca, kmorris@uwaterloo. ca 


Abstract 

The Landau-Lifshitz equation describes the dynamics of magnetization inside a ferromagnet. This equation is nonlinear and 
has an infinite number of stable equilibria. It is desirable to control the system from one equilibrium to another. A control 
that moves the system from an arbitrary initial state, including an equilibrium point, to a specified equilibrium is presented. 
It is proven that the second point is an asymptotically stable equilibrium of the controlled system. The results are illustrated 
with some simulations. 
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1 Introduction 

The Landau-Lifshitz equation describes the magnetic 
behaviour within ferromagnetic structures. This equa¬ 
tion was originally developed to model the behaviour of 
domain walls, which separate magnetic regions within a 
ferromagnet [1]. Ferromagnets are often found in mem¬ 
ory storage devices such as hard disks, credit cards or 
tape recordings. Each set of data stored in a memory 
device is uniquely assigned to a specific stable magnetic 
state of the ferromagnet, and hence it is desirable to con¬ 
trol magnetization between different stable equilibria. 
This is difficult due to the presence of hysteresis in the 
Landau-Lifshitz equation. Due to the presence of mul¬ 
tiple equilibria, a particular control can lead to different 
magnetizations. The particular path depends on the ini¬ 
tial state of the system and looping in the input-output 
map is typical [2,3]. 

There is now an extensive body of results on control 
and stabilization of linear partial differential equations 
(PDE’s); see for instance the books [4,5,6,7] and the re¬ 
view paper [8]. There are much fewer results on control 
and stabilization of nonlinear PDE’s and the Landau- 
Lifshitz equation is particularly problematic. Experi¬ 
ments and numerical simulations demonstrating the con¬ 
trol of domains walls in a nanowire are presented in 
[9,10]. In [11], the Landau-Lifshitz equation is linearized 
and shown to have an unstable equilibrium. One control 
objective is to stabilize this equilibrium with the control 
as the average of the magnetization in one direction and 
zero in the other two directions. In [12,13], solutions to 
the Landau-Lifshitz equation are shown to be arbitrar¬ 


ily close to domain walls given a constant control. 

Stability results are often based on linearization 
[14,15,16,17]. However, as is well known, results with a 
linearized model can only predict local stability. Fur¬ 
thermore, for PDFs models, stability of the linearization 
does not necessarily predict even local stability of the 
original model [18]. The nonlinearity in the Landau- 
Lifshitz equation is discontinuous, that is, the equation 
is not quasi-linear, and there are no results that can 
be used to conclude that local stability of the nonlin¬ 
ear equation follows from the nonlinear equation [19]. 
Furthermore, in many applications the goal is to move 
the system from one equilibrium to another equilibrium 
and linearization is not always useful in this context. 

In the next section, the uncontrolled Landau-Lifshitz 
equation is described. The main result on stabilization 
of the Landau-Lifshitz equation is in Section 3. In Sec¬ 
tion 4, simulations illustrating the results are presented. 


2 Landau-Lifshitz Equation 

Consider the magnetization 

ni(x,t) = (toi(x, 1), TO 2 (x, 1), TO 3 (x, 1)), 

at position x € [0, L] and time t > 0 in a long thin ferro¬ 
magnetic material of length L > 0. If only the exchange 
energy term is considered, the magnetization is modelled 
by the one-dimensional (uncontrolled) Landau-Lifshitz 
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equation [20],[21, Chapter 6 ] 

^ = m X - i/m X (m X (la) 

ot 

m(a;, 0 ) = mo (x) (lb) 

where x denotes the cross product and > 0 is the 
damping parameter, which depends on the type of fer- 
romagnet. The Landau-Lifshitz equation sometimes in¬ 
cludes a parameter called the gyromagnetic ratio multi¬ 
plying m X m.xx ^ where m^x means the magnetization is 
differentiated with respect to x twice. This has been set 
to 1 for simplicity. For more on the damping parameter 
and gryomagnetic ratio, see [ 22 ]. 

The Landau-Lifshitz equation is a coupled set of three 
nonlinear PDFs. It is assumed that there is no magnetic 
flux at the boundaries and so Neumann boundary con¬ 
ditions are appropriate: 

ma;( 0 , t) = irixiL, t) =0 (Ic) 


Ferromagnets are magnetized to saturation [26, Sec¬ 
tion 4.1]; that is ||mo(x )||2 = Mg where || • [[2 is the 
Euclidean norm and Mg is the magnetization satura¬ 
tion. In much of the literature. Mg is set to 1 ; see for 
example, [21, Section 6.3.1], [23,24,27]. This convention 
is used here. Physically, this means that at each point, 
X, the magnitude of mo(x) equals the magnetization 
saturation. The initial condition mo(x) is furthermore 
assumed to be real-valued. The assumption on the ini¬ 
tial magnetization is satisfied by the magnetization at 
all time. 

The set of equilibrium points of (1) is [21, Theorem 6.1.1] 
E ={a = (tti, 02 , 03 ) : Oi, 02 , 03 constants and a^'^a = 1}. 

( 5 ) 

Theorem 4. [25, Theorem 4.11] The equilibrium set in 
(5) is asymptotically stable in the C 2 ~''^orm. 

3 Controller Design 


where m^; means the magnetization is differentiated with A control, u(f), is introduced into the Landau-Lifshitz 
respect to x once. equation (la) as follows 


Existence and uniqueness of solutions to (1) with differ¬ 
ent degrees of regularity has been shown [23,24]. 

Theorem 1. [21, Lemma 6.3.1] //||mo(x )||2 = 1, the 
solution, m, to (la) satisfies 


m(x,t)||2 = 1. 


( 2 ) 


9m 

—— = m X nxixx 

dt 

m(x, 0 ) = mo(x). 


I'm X (m X mxx) + u(t) ( 6 ) 


As for the uncontrolled system, the boundary condi¬ 
tions are ma;(0,t) = mx{L,t) = 0. Equation (6) is 
the Landau-Lifshitz equation with an external magnetic 
field u. 


The following statement is a more restrictive version of 
the theorem stated in [23]. 

Theorem 2 . [23, Thm. 1.3,1.4]- If rno € H 2 { 0 ,L), 
mo,£!;( 0 ) = mo,a:(T) = 0 and Umo]] = 1 , then there exists 
a time T* > 0 and a unique solution m of (1) such that for 
all T <T*, me C([0, T; 772 ( 0 , L)) n £ 2 ( 0 , 7; 7 / 3 ( 0 ,7)). 

With more general initial conditions, solutions to (1) 
are defined on £2 = 7i2([0, 7]; M^) with the usual inner- 
product and norm. The notation || • ||£3 is used for the 
norm. Define the operator 

/(m) = m X mxx -vmx{mx mxx), (3) 

and its domain 

D = {meCl:mxe L\, mxx G 

m,j,(0) = m,j.(7) = 0}. (4) 

Theorems. [25, Theorem 4.7] The operator f (m) with 
domain D generates a nonlinear contraction semigroup 
on £ 3 . 


The goal is to choose a control so that the system 
governed by the Landau-Lifshitz equation moves from 
an arbitrary initial condition, possibly an equilibrium 
point, to a specified equilibrium point r. The control 
needs to be chosen so that r becomes a stable equilib¬ 
rium point of the controlled system. It can be shown 
that zero is an eigenvalue of the linearized uncontrolled 
Landau-Lifshitz equation [25, Chapter 4.3.2]. For finite¬ 
dimensional linear systems, simple proportional control 
of a system with a zero eigenvalue yields asymptotic 
tracking of a specified state and this motivates choosing 
the control 

u = 7(r — m) (7) 

where r € 71 is an equilibrium point of the uncontrolled 
equation ( 1 ) and fc is a positive constant control param¬ 
eter. 

The following theorem shows that for any initial condi¬ 
tion mo the solution to ( 6 ) with control u{t) = k{r — m) 
satisfies ||m {■,t)\\cl < 1 - 

Theorem 5. For any r G E define 

Bm = k{r — m). ( 8 ) 
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If k > 0, the nonlinear operator f + B with domain 
D, where f and D are defined in (3), (4) respectively, 
generates a nonlinear contraction semigroup on L\. 


PROOF. It will be shown that (i) f + B is dissipative 
and (ii) the range of I — a{f + B) for all a > 0 is £ 2 - 

For any m.,y G D, 

(/(m) + Bm - (/(y) + By), m - y)^ 

= (/(m) - /(y), m - y)c3 + {Bm ~ By, m- y)c3. 

Since / generates a nonlinear contraction semigroup 
(Theorem 3), / is dissipative [28, Proposition 2.98] and 
hence 

(./■(m)-/(y),m-y)£3 < 0. 

It follows that 

(/(m) + Bm - (/(y) + By),m - y) ^ 

< {Bm-By,m-y)ci 

= {-km + ky,m-y)^3 

= -k{m-y,m-y)ci 

< 0 


The results in Lemma 6 and 7 are a consequence of the 
product rule. 

Lemma 6. For m G C\, the derivative o/g = m x m^, 
is g:r = m X m^x ■ 

Lemma 7 . For m G C2, the derivative of f = 
(m X ma,)"'" (m x m^,) is fx =2 (m x m^,)"^ (m x mxx) ■ 

Lemma 8. Form G £3 satisfying (Ic), 

/ (m — r)"'"(m x mxx)dx = 0. 

Jo 


PROOF. Integrating by parts, and applying Lemma 6 
and the boundary conditions (Ic) implies that 

pL pL 

/ (m — r)"'"(m x mxx)dx = — mj(m x mx)dx. 

Jo Jo 

From properties of cross products, mj(m x m^,) = 
m"^(m 2 , X ma;) = 0, and hence the integral is zero. □ 

Lemma 9. Form G £2 satisfying (Ic), 


and hence / + F is dissipative. 


m X ma ;||£3 < 4L^||m x m 


xx\\c% 


Since / generates a nonlinear contraction semigroup 
(Theorem 3), then it is m-dissipative and hence, 
ran(I — af) = C\ for any d > 0 [29, Lemma 2.1]. This 
means that for any y^ G C\ there exists m G F such 
that m — d/(m) = y^- Choose any yi G £ 2 , a > 0 and 
define 

yi akv 

y 2 = 


1 + afc 1 + abk 


and 


1 + ak 

There exists m G F such that 
/(m) = y 2 = 


yi 


Offer 


1 + ofe 


I + ofe 1 + abk 


PROOF. Integrating by parts, using Lemma 7 and the 
boundary conditions (ic) leads to 


m X m 


■x\\c% 


T 

/ 2(mxma,) {mxmxx)xdx. 

Jo 


It follows from Young’s inequality that 

cL 


1 

m X ma;||^3 ^ 7 : X X dx 

^ 2 Jo 

fL 

+ 2 (m X mxx) (m x ^xx) x^dx. 

Jo 


Solving for yi leads to 

yi = m — a(fe(r — m) + /(m)). 

Thus, for any yi G £3, there exists m G F such that 
yi = (I—a(F+/))mandhenceran(I — a{B + /)) = £3 
for some o > 0. It follows that the range is £| for all 
o > 0 [29, Lemma 2.1]. 

Thus, since F+/ is dissipative and the range of (I—o(F+ 
/)) is £ 2 , then B + f generates a nonlinear contraction 
semigroup [28, Proposition 2.114]. □ 


Since x G [0, F], 

||m X m,j;||^3 < ^||m X mx\\\3^ + 2L^\\m x mxxW^i^- 

Rearranging gives the desired inequality. □ 

Theorem 10. Let r be an eguilibrium point of (6). For 
any constant k such that k > 8vL‘^, r is a globally expo¬ 
nentially stable equilibrium point of (6) in the Hf-norm. 
That is, for any initial condition on Hf, m decreases ex¬ 
ponentially in the Hi-norm to r. 
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PROOF. If the initial condition mo{x) G Hf, then so¬ 
lutions to ( 6 ) are in Hf [24]. The Lyapunov candidate 

F(m) = ^ ||m-r||^3 + ^ ||m,||^3 

is clearly positive definite for all m G D. Furthermore, 
V (m) = 0 only when m = r. Taking the derivative of V 
along trajectories of the controlled equation ( 6 ) 


— = / (m — -f / 

dt Jo Jo 

pL p L 

= / (m — v)^Yadx — / rn^^uidx 

Jo Jo 

where the dot notation means differentiation with re¬ 
spect to t. Substituting in (6) to eliminate rh, 

^ = J (m - r)'^ (m X dx 

— ^ (m — (m x (m x nia^a^)) dx 

Jo 

-\-k (m — r)"’"(r — m)(ix 
^0 

- / (m X m.^x)dx 

Jo 

+ (m X (m x dx 

Jo 

m^x{^-m)dx. 

Jo 

From Lemma 8, the first integral is zero. Furthermore, 
from properties of cross products, 

Tcn^^ (m X 013,3,) = (103,3, x 013,3,) = 0 , 


and hence 


/ (01 X 013 , 3 ,) dx = 0. 

Jo 


It follows that 


^ = - V J (01 - r)'’" (m X (m X 013 , 3 ,)) dx 

- fcljm - r ||£3 - i^ljoi X 013,3,11^3 
•^2 ^2 

-k m.l^{r-m)dx. 

Jo 

Applying integration by parts to the last integral, 
[ m)da; = ||oi3,||£3. 

Jo 


The first integral in ^ can be rewritten 
—v / (m — r)"*" (m x (01 x 013,3,)) dx 

Jo 

= — V j ((m — r) X m) (m x 103,3,) dx 

Jo 

X -P 

= 1 ^ / (r X 01) (m X 013,3,) dx 

Jo 

Applying integration by parts to this integral, using 
Lemma 6 and boundary conditions (Ic) leads to 

X P 

v (r X m) (m x 013,3,) dx = —v{y x 0I3,, m x oi3,)£3. 

Jo "" 

Then from Cauchy-Schwarz, Lemma 9, Young’s Inequal¬ 
ity and Lemma 7, 

X P 

V / (r X m) (m X 013,3,) dx < jzjjr X 013,11 £3II01 X 013,11 £3 

Jo 

< 4 j^L^||r X oi 3 ,||£ 3 ||oi X oi 3 , 3,||£3 

< 8 :^L'‘||oi 3,||£3 -f ^lloi X Ol 3 , 3 ,||£ 3 . 

Substituting this inequality and (10) into (9) leads to 

- ^Ilm X 013,3, 1 1£3 - A:||oi - r||£3 
< — (k — ||oi3,||£3 — A:||oi — r||£3 

<-{k- (^||oi 3,||£3 + llm- r||£ 3 ^ 

= -2 (fc - V. 

Integrating with respect to time and noting that r does 
not depend on x 


< (||m,(x, 0 ) 11^3 + ||m(x, 0 ) -r||23) 

= (11 (m(x,0) - r), 11^3 + ||m(x,0) - rH^^) . 


Therefore, 


<e- 2 (fc- 8 -i")*||m(x, 0 )-: 


and since k > 8vL^, r is an exponentially stable equilib¬ 
rium point of ( 6 ). □ 


A natural question is whether r is exponentially stable 
in the Analysis of the linear Landau-Lifshitz 


4 



equation provides insight to this question. The linearized 
controlled Landau-Lifshitz equation is 

— = + ax + k{r - z) , z(0 ) = Zq (11 ) 

with the same boundary conditions Z 2 ;( 0 ) = Zx{L) = 0 . 
Since the uncontrolled linear Landau-Lifshitz equation 
generates a linear semigroup and /c (r — z) is a bounded 
linear (affine) operator, then the operator in (11) gener¬ 
ates a semigroup [5, Theorem 3.2.1]. Substituting z = r 
into ( 11 ) leads to dz/dt = Q and hence r is a stable equi¬ 
librium point of ( 11 ). 

Theorem 11 . Let r G E. For any positive constant k, 
r is an exponentially stable equilibrium of the linearized 
system (11) in C^-norm. 


PROOF. For z G D{A), where D{A) = H as in equa¬ 
tion (4), consider the Lyapunov candidate 

^(z) = ^I|z-i-||£3. 

It is clear that F > 0 for all z G D{A) and furthermore, 
V {z) =0 only when z = r. Therefore, V{z) >0 for all 
z G L'(^)\{r}. 

Taking the derivative of V (z) implies 
^ = I (z - r)^z<ix. 


Solving yields 




< e 


— 2kt I 


zo - r 


For fc > 0 the equilibrium point, r, of (11) is locally 
exponentially stable. This is true for any initial condition 
and hence global stability is obtained. □ 


Theorem 11 suggests that the equilibrium point in the 
controlled nonlinear Landau-Lifshitz equation ( 6 ) is 
exponentially stable in the T^-norm. However, since 
the nonlinearity in the Landau-Lifshitz equation is un¬ 
bounded, stability of the linear equation does not neces¬ 
sarily reflect stability of the original nonlinear equation; 
see [18,19]. 


In equation ( 6 ), the control is affine. However, in current 
applications, the control enters as an applied magnetic 
field [11,12,13,14,15]. More precisely. 


(9m 


=m X {nixx + u)- izmx {mx {nixx + u)) 

=m X nixx — x (m x nixx) 

-I- m X u — izm X (m x u) (12) 


where u = fc(r—m) as before and r G i? is an equilibrium 
point of (12). Equation (12) is the Landau-Lifshitz equa¬ 
tion with a nonlinear control. Its existence and unique¬ 
ness results can be found in [23, Thm. 1.1,1.2] and is 
similar to Theorem 2. 


As for the uncontrolled equation, since 


Substituting in (11) yields 


— = v j (z - Y)'^Zxxdx + (r - z)'^ (a x z^x) dx 
Jo Jo 

+ k {z — r)'^{r — z)dx. 

Jo 

By Lemma 8 , the middle term is zero. Using integration 
by parts, the first term becomes 

T 

—ly / zfzxdx. 

Jo 

It follows that 


1911m(x,t)ll2 rpclm 

-=m - 

2 dt dt 

=m'^(m X m.xx — x (m x nixx) 

-I- m X u — I'm X (m x u)) = 0, 

this implies ||mjj 2 = c, where c is a constant. The con¬ 
vention is to take c = 1. It follows that any equilibrium 
point is trivially stable in the >C 2 -norm. 

Theorem 12. For any r G E and any positive constant 
k, Y is a locally stable equilibrium point of (12) in the 
Hi-norm. That is, for any initial condition mo{x) G D, 
where D is defined in (4), the Fli-norm of the error m — Y 
does not increase. 


dt 


-i'\\zx\\'^^3 - k\\z 



and since v > 0, 


dt 


< —k\\z — rjj ^3 = —2kV. 


PROOF. Let B{y,p) = {m G £2 ■ 11^^“ ''llcl < p} C 
D for some constant 0 < p < 2. Note that since p < 2, 
then —r ^ B(y,p). For any m G B{y,p), consider the 
iLi-norm of the error 

U(m) = fcllm- rll^i + WnixWls . 
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Taking the derivative of V, 


PROOF. For an admissible v with | |v — r| I 2 < 2, con¬ 
sider the Lyapunov candidate 


dV 

dt 



k{m—Y)^mdx+ / 

Jo 

fc(m — Y)^mdx — / mj^rhda; 

Jo 

(fc(m — r)"'"m — m^^m) dx. 


(13) 


y(v) 



It is clear that F > 0 for all v e D{A) and furthermore, 
V (v) = 0 only when v = r. Therefore, V (v) > 0 for all 
V e I?(A)\{r}. 


Let h = m — r, then the integrand becomes 

fch'^m — mj^rh (14) 

and equation (12) becomes 

m = m X (nixx — Ah) — lym x (m x (nHxx — Ah)). 

It follows that 

h'^rh =h'^ (m x nixx) — v (m x ma,^;)'^ (h x m) 

— i/A||m X h ||2 (15) 

and 


Taking the derivative of V (v) implies 
/ (v — r)"'"vdx 

Jo 

and substituting in (17) leads to 
dV 


dV 

dt 


dt 


plx p L/ 

V {■v-r)'^Vxxdx + iy {v - r)'^(r x ■Vxx)da 
Jo Jo 

+ k (v — r)'^(v X r)da; 

Jo 

— kiy (v — r)"'’ (r X (v X r)) dx. 

Jo 


mJajih = — knixx (m x h) + v\\u\ x ma;a ;||2 

+ vk (m X h)"^ {too-xx x m). (16) 


The second integral is zero by Lemma 8, and applying 
integration by parts the first term becomes 


Substituting (15) and (16) into equation (14) leads to 
fch'^m — mja;ih = 2 vk (m x m.xx)^ (m x h) 

— ^A^||m X h||2 — i^||m X ma;a;||2 

= — ^||m X nixx — Am x h||2 

Substituting this expression into equation (13) leads to 

dlX \ 11 2 l■\ 

— = - v\\m X [nixx + u)||£3 < 0. 

Thus, the iJi-norm of the error does not increase. □ 


For any equilibrium point r € i? of (12) and m G D, let 
m = r + V where v is any admissible perturbation; that 
is V G H and ||r + v ||2 = 1. The linearization of (12) at 
r is 


9v 


Wxx 


v(0 ) = Vq. 


r X Vxx + Av X r — ki'r x (v x r) 


(17) 


—ly / vjvxdx. 


It follows that 


dV 

dt 


pL pL 

= — vi if'^'Vxdx + k / (v — r)"'"(v X r)da; 

Jo Jo 

— kv j (v — r)"^ (r X (v X r)) dx 

Jo 

3 + A / h'^(h X r)dx 

^ Jo 

— kv j (r X (h X r)) dx 

Jo 


= - V; 


where h = v — r. The middle integral is zero since 
h"’’(h X r) = r"^(h X h) = 0, and the last integral can be 
simplified using the fact that 

(r X (h X r)) = (h x r)"^ (h x r). 


Therefore, 


dV 

dt 


= -J^lWx 


2 

^2 


ki 2 \\h X rll^s. 


Theorem 13. Letr G E be an equilibrium point of (17). 
For any positive constant k, r is a locally asymptotically 
stable equilibrium of (17) in the C^-norm. 


6 



For k > 0, 

^ =-v (||v ^||£3 + fc||h X < 0 

and furthermore, dV/dt = 0 if and only if v^, = 0 and 
hxr = vxr = 0. This is true only if v = ar where 
a is any scalar. Since v must be a constant satisfying 
V X r = 0, then v = ar for some constant a. But since 
||r + v ||2 = 1 , then a = 0 . 

It follows that r is a locally asymptotically stable equi¬ 
librium point of (17). □ 


4 Example 

Simulations illustrating the stabilization of the Landau- 
Lifshitz equation were done using a Galerkin approxi¬ 
mation with 12 linear spline elements. For the following 
simulations, the parameters are v = 0.02 and L = 1 
with initial condition mo (a;) = (sin(27rx), cos(27rx), 0). 
Figure 1 illustrates that the solution to the uncontrolled 
Landau-Lifshitz equation settles to ro = (0, —0.6,0). 

Stabilization of the Landau-Lifshitz equation with affine 
control ( 6 ) is illustrated in Figure 2 with the second 
equilibrium point chosen to be ri = (—^,0, ^). The 
control parameter is fc = 0.5. Figure 3 depicts applying 
the control twice in succession, forcing the system from 
the equilibrium Tq to r 2 = ( 1 , 0 , 0 ) and then to a new 
equilibrium ra = (0, 0,1). In each case, the state of the 
controlled system converges to the specified point as 
predicted by the analysis. 

Stabilization of the Landau-Lifshitz equation with non¬ 
linear control (12) is illustrated in Figure 4 with the sec¬ 
ond equilibrium point chosen to be ri. The control pa¬ 
rameter is k = 10. It is clear from the figure that the 
system converges to the specified equilibrium point, ri. 
The control can also be applied after the dynamics have 
settled to Tq as shown in Figure 5. In Figure 6 , the dy¬ 
namics settle (without a control) to rg, then the control 
is applied in succession twice, which forces the system 
from rg to ri, and then finally to r 4 = (0,1, 0). The rate 
of convergence is slower and the value of k needed is 
larger than with an affine control. 

5 Conclusion 

The Landau-Lifshitz equation is a nonlinear system of 
PDFs with multiple equilibrium points. The fact that 
it is not quasi-linear means that the linearization is not 
guaranteed to predict stability of the non-linear equa¬ 
tion [19]. Furthermore, since the objective of the control 
is to steer between equilibrium points, a linearized anal¬ 
ysis, which only yields local results, would not predict 


stability of the controlled system. However, the presence 
of a 0 eigenvalue in the linearized equation suggested 
that simple feedback proportional control could be used 
to steer the system to an arbitrary equilibrium point. 

The controlled system with an affine control term was 
shown to be well-posed and also globally asymptotically 
stable. Furthermore, it is exponentially stable in the Hi- 
norm. Analysis of the linearized Landau-Lifshitz equa¬ 
tion, which shows that the equilibrium is exponentially 
stable, suggests the original nonlinear Landau-Lifshitz 
equation is locally exponentially stable in the L 2 -norm. 
Future research aims to establish this. 

In applications, the control enters through an applied 
field and the control enters nonlinearly. It was shown 
that the Landau-Lifshitz equation with a nonlinear con¬ 
trol has a stable equilibrium point and the linearization 
has an asymptotically stable equilibrium point. Simu¬ 
lations indicate that proportional control also stabilize 
the fully nonlinear model. This suggests the nonlinear 
equation has an asymptotically stable equilibrium point. 
Proving this remains an open research problem. 
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Fig. 1. Magnetization in the uncontrolled Landau-Lifshitz 
equation moves from initial condition mo (a:), to the equilib¬ 
rium ro = (0, —0.6,0). 


i> 1 



Fig. 2. With a proportional control {k = 0.5), magnetization 
in the Landau-Lifshitz equation with a linear control moves 
from the initial condition mo (a:) to the specified equilibrium 







































Fig. 3. Steering magnetization between specified equilibria 
with a linear control. The uncontrolled magnetization moves 
from initial condition mo to ro = (0,—0.6,0). Proportional 
control {k = 0.5) with two successive values of r first forces 
the magnetization to r 2 = (1,0, 0) and then to ra = (0,0,1). 


Fig. 5. Magnetization in the Landau-Lifshitz equation 
moves from the initial condition mo (a;) to the equilibrium 
To = (0, —0.6,0) without a control. The control, u with 
fe = 10 is then applied to the equation nonlinearly and steers 
the dynamics to the specified equilibrium ri = (“^’0’^)- 




Fig. 4. Magnetization in the Landau-Lifshitz equation with 
nonlinear control moves from the initial condition mo (a:) 
to the specified equilibrium ri = 1^) ''''ith control 

parameter k — 10. 





Fig. 6. Magnetization in the Landau-Lifshitz equation 
moves from the initial condition mo (a;) to the equilibrium 
To = (0, —0.6,0) without a control. The control u with 
fe = 10 is then applied to the equation nonlinearly and steers 
the dynamics to the specified equilibrium ri = O’ 

and then to another equilibrium, r 4 = (0,1, 0). 
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